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Abstract 

We give a closed form for the correlation functions of ensembles of a class of asym- 
metric real matrices in terms of the Pfafnan of an antisymmetric matrix formed from a 
2x2 matrix kernel associated to the ensemble. We apply this result to the real Ginibre 
ensemble and compute the bulk and edge scaling limits of its correlation functions as 
the size of the matrices becomes large. 
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1 Introduction 

The principal subject of this paper is the Ginibre ensemble of real random matrices - 
square real valued matrices with i.i.d. normal entries. Random matrix models have been 
very successful in describing various physical phenomena. (See e.g. [16J and references 
therein). Physical applications of the Ginibre ensembles are described in [13], [17] and 
|12j . Mathematically, the ability to analyze spectra of random matrices is largely based 
on determinental or Pfaffian formulas for spectral correlations that have been derived for a 
variety of models. The real Ginibre ensemble is one of the few models for which such formulas 
remained unavailable for over 40 years (the model was introduced by Ginibre in 1965 [M]). 
The goal of this paper is to prove Pfaffian formulas for all the correlation functions of the 
real Ginibre ensemble and to evaluate their bulk and edge scaling limits. 

The algebraic techniques that we use were developed in a recent work of the second 
author [21]. Starting from that paper, Forrester and Nagao [12] have independently ob- 
tained similar Pfaffian formulas for correlations and constructed certain skew-orthogonal 
polynomials necessary for the asymptotic analysis. We take the next step and obtain the 
asymptotics of the correlation functions using these polynomials. 

In the algebraic part of this paper we consider a general class of probability measures that 
includes the real Ginibre ensembles and ensembles arising in the study of Mahler measure 
of polynomials; the latter are of interest in number theory, see [22]. We show that the 
correlation functions for an ensemble from this class can be expressed as the Pfaffian of 
a block matrix whose entries are expressed in terms of a 2 x 2 matrix kernel associated 
to the ensemble. We find much inspiration from Tracy and Widom's paper on correlation 
and cluster functions of Hcrmitian and related ensembles [26] . However, instead of using 
properties of the Fredholm determinant to calculate the correlation functions via the cluster 
functions, we use the notion of the Fredholm Pfaffian to determine the correlation functions 
directly. A Pfaffian analog of the Cauchy-Binet formula introduced by Rains [20] lies at 
the heart of our proof. For completeness we will include Rains' proof here. In place of the 
identities of de Bruijn [7] used by Tracy and Widom we will use an identity of the second 
author [21j to compute the correlation functions for ensembles of asymmetric matrices. 
Rains' Cauchy-Binet formula has applicability in a wider context than just the determination 
of the correlation function for ensembles of asymmetric matrices, and we will use it to give 
a simplified proof of the correlation function of Hermitian ensembles of random matrices 
when (3=1 and (3 = 4. 

We then apply our theory to the real Ginibre ensemble. It is known (see [15] [2] [9] 
[5]) that the density functions of real and complex eigenvalues is approximately constant on 
(— >/~N, yN) and the disc of radius y/~N respectively (here N is the size of the matrices). Wc 
study the local correlations of eigenvalues in four different regions: near a real point u^f~N 
with — 1 < u < 1 (real bulk), near ±\/N (real edge), near a (non-real) complex point u\/~N 
with |u| < y/~N (complex bulk), and near u^/N with |u| = y/N, Im(w) ^ (complex edge). 
Four different limit processes arise, and we compute their correlation functions explicitly. 
The complex bulk and edge limits turn out to be the same as in the case of the much simpler 
complex Ginibre ensembles. The correlations of real eigenvalues in the real bulk region were 
obtained in |12] and the density functions were computed earlier using different techniques 
[8] [9]. All other results appear to be new. 

The paper is organized as follows. Section [2] introduces necessary notation. In Section[3j 
we introduce a class of ensembles relevant for our study. In Section [4] we show that the 
real Ginibre ensemble and the ensemble related to Mahler measure fall into this class. In 
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Section [51 we introduce the correlation functions. In Section [6] we construct the correlation 
kernel, and in Section [7] we state how the correlation functions are expressed through that 
kernel. Section[5]contains statements of asymptotic results as well as the limiting correlation 
kernels. Section [H] contains the proofs. Appendix [A] shows how to apply Rains' Cauchy- 
Binet formula to the (3 = 1, 4 Hermitian random matrix ensembles. Appendix [Blcontains the 
proof of Rains' formula. In Appendix [Cl we compute the bulk and edge limits of the complex 
Ginibre ensemble. (We were unable to locate this material in the literature). Appendix [Dl 
provides plots of the first and second correlation functions in various limit regimes. 

Acknowledgments. We thank Percy Dcift for helpful early discussions and for introducing 
the authors to each other. The second author would also like to thank Brian Rider for many 
helpful discussions regarding the Ginibre real ensemble. We are grateful to Eric Rains for 
allowing us to include his proof of the the Pfaffian Cauchy-Binet formula (Appendix [B| . 
Finally we would like to thank Peter Forrester for keeping us updated about his work. The 
first named author (A.B.) was partially supported by the NSF grant DMS-0707163. 

2 Point Processes on the Space of Eigenvalues of R iVx7V 

We start rather generally since the results in this manuscript can be used to describe not 
only the statistics of eigenvalues of ensembles of real matrices but also the statistics of 
roots of certain ensembles of real polynomials. We begin with random point processes on 
two-component systems. 

Let X be the set of finite multisets of the closed upper half plane H C C. An element 
£ G X is called a configuration, and X is called the configuration space of H. Given a 
Borel set A of H, we define the function Na ■ X — ► Z-° by specifying that Xa(£) be the 
cardinality (as a multiset) of (£ PI A). We define £ to be the sigma-algebra on X generated 
by {Na : A C H Borel}. A probability measure CP defined on £ is called a random point 
process on X. 

For each pair of non-negative integers (L,M) we define Xl^m to be the subset of X 
consisting of those configurations £ which consist of exactly L real points and M points in 
the open upper half plane H. That is, 

Xl,m ■■= {i G X : Nr(0 = L and N ff (0 = Af}. 
Clearly, Xl,m is measurable and X can be written as the disjoint union 

A = I^J Xl^m- 

L>0,M>0 

Given a point process CP on X, we may define the measure Tl.m on X by CPj^jw^-B) := 
7{B n Xl,m) f° r each B G E . The measure Tl.m induces a measure on Xl : m (and we will 
also use the symbol Vl,m for this measure). 

Given a matrix Y G M. NxN there must be a pair of non-negative integers (L,M) with 
L + 2M = N such that, counting multiplicities, Y has L real eigenvalues and M non-real 
complex conjugate pairs of eigenvalues. By representing each pair of complex conjugate 
eigenvalues by its representative in H, we may identify all possible multisets of eigenvalues 
of matrices in ~R NxN with the disjoint union 

X N := [J A L , M . (2.1) 

(L,M) 
L + 2M = N 
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Similarly, we may identify all possible multisets of roots of degree N real polynomials with 
this disjoint union. Thus, when studying the statistics of eigenvalues of ensembles of asym- 
metric real matrices (respectively of the roots of degree N real polynomials) we may restrict 
ourselves to random point processes on H which are supported on the disjoint union given 
in (|2.ip . That is, the eigenvalue statistics of an ensemble of real matrices is determined 
by a set of finite measures T'l.m on Xl,m for each pair of non-negative integers (L,M) 
with L + 2M = N. In this situation we will say that 7 is a random point process on Xn 
associated to the family of finite measures {!Pl,m : L + 2M = N}. 

From here forward we will assume that L, M and N are non- negative integers such that 
L + 2M = N, and a sum indexed over {L, M) will be taken to be a sum over all pairs 
satisfying this condition. Moreover, when we refer to a point process we will always mean 
a point process on X^. 

3 Point Processes on Associated to Weights 

In this section we will introduce an important class of point processes associated to Borcl 
measures on C. We will be particularly interested in measures which are a sum of two 
mutually singular measures: one of which is absolutely continuous with respect to Lebesgue 
measure on R and one which is absolutely continuous with respect to Lebesgue measure on 
C. The corresponding densities with respect to Lebesgue measure will allow us to construct 
a weight function which uniquely determines the associated point process on Xn- Point 
processes of this sort arise in the study of asymmetric random matrices and the range of 
multiplicative functions on polynomials with real coefficients. 

It will be useful to distinguish non-real complex numbers and we set C* = C \ R. 

We start rather generally by constructing measures on Xl,m from measures on R L x C^ 1 . 
The benefit in doing this is that it allows us to express important quantities associated to 
point processes (averages, correlation functions, etc.) as rather pedestrian integrals over 
R L x C^ 1 . To each (a,/3) € R L x we associate a configuration in Xl,m given by 

{c*,/3} :={«!,... ,£*£,/?!, ... 3m}, where [3 m = {/3 m , f3 m } n H. 

A given configuration £ S Xl.m may correspond to several vectors in R L x and we will 
call {(a, pi) : {a,/3} = £} the set of configuration vectors of £. 

A function F on Xl m induces a function on R L x Cf specified by (a,0) i-> F{a,f3}, 
and given a measure vl,m on R L x there exists a unique measure 7l,m on Xl,m specified 
by demanding that 

f F(Z)dP L , M (t)= n ,L M / F{ a ,[3}du L>M (a,f3), (3.1) 

JX L>M hUVlll Js. L xCy 

for every X-measurable function F on X. The normalization constant LIM\2 M arises since 
a generic element £ € Xl,m corresponds to L\M\2 M configuration vectors. By specifying 
measures vl.m on R L x for all pairs (L, M) and normalizing so that the total measure 
of Ajv is 1, we define a point process on Ajy. 

A very important class of point processes arises when we demand that the various vl.m 
are all related to a single measure on C. Given a measure v\ on R and a measure v-i on 
C*, we set v to be the measure v\ + v<i on C. We will write vl for the product measure 
of v\ on R L and V2M will be the product measure of on < C M . By combining vl and 
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Vim with a certain Vandermonde determinant we will arrive at the desired measures on 
R L x C* f . Given a vector 7 e C N we define V(*y) to be the N x N matrix whose n, n' entry 
is given by 7^ _1 . We will denote the determinant of V{~f) by A(7), and define the function 
A : R L x C M -> C by 

A(a,/3) := A(ai, . . . , az,,/3i,/3i, . . . ,Pm,(3m). 
Using these definitions we set isl.m to be the measure on R L x given by 

f" L ' M («,/3)=2^|A(a,/3)l, (3.2) 

and we will write CP£ M for the measure on Xl,m given specified by vl,m as in (|3.f p . If 
CP£ M is finite for each pair (L, M), then we set 

(L,M) (L.M) 

We will call CP 1 ' the point process on Xp? associated to the weight measure v and Z v will be 
referred to as the partition function of y v . 

We set Ai and A2 to be Lebesgue measure on R and C respectively and let A := Ai + A 2 . 
If there exists a function w : C — > [0, 00) such that v — wX (by which we mean dv jdX = w), 
then we will call w the weight function of CP". In this situation we set M := CP£ M and 
define J 11 " and Z w analogously. If P" has weight function w then 

, , h M 

dv L ,M 



d(X L x A 2 m) 



II H3m)}|A(a,/3)|, (3.3) 

£=1 m=l ' 



and we define CIl,m '■ K L x C M to be the function given on the right hand side of f|3 .3[) . 
The collection {Ql,m ■ L + 2M — N} plays the role of the joint eigenvalue probability 
density function, and we will call Ol,m the L, M- partial joint eigenvalue probability density 
function of CP™. 



4 Examples of Point Processes Associated to Weights 

Point processes on Xn associated to weights arise in a variety of contexts. It is often the 
case that the weight function w is invariant under complex conjugation. In this situation, 
there necessarily exists some function p : C — > [0, 00) such that 

r p( 7 ) if 7 GR; 
Ml) : = I 

[ p{i)p{l) if 7 e C*. 

When it exists, we will cal p the root or eigenvalue function of CP™ (depending on whether 
CP™ models the roots of random polynomials or eigenvalues of random matrices). The 
root /eigenvalue function is often a more natural descriptor of CP™ than w. In this situa- 
tion we will write CP^ and for 7 W and Z w . 
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4.1 The Real Ginibre Ensemble 

In 1965 J. Ginibre introduced three ensembles of random matrices whose entries were re- 
spectively chosen with Gaussian density from K, C and Hamilton's quaternions 14]. The 
real Ginibre ensemble is given by M. NxN together with the probability measure ij given by 

d V (Y) = (2ny N2 ^e-^ T ^dX NxN (Y) 7 

where Xnxn is Lebesgue measure on M. NxN . This ensemble has since been named Gi- 
nOE due to certain similarities with the Gaussian Orthogonal Ensemble. Among Ginibre's 
original goals was to produce a formula for the partial joint eigenvalue probability density 
functions of GinOE. He was only able to do this for the subset of matrices with all real 
eigenvalues. In the 1990s Lehmann and Sommers [TH] and later Edelman [5] proved that 
the L, M-partial eigenvalue probability density function of GinOE is given by £Il,m for the 
eigenvalue function 

Gin(7) := exp(— 7 /2) V /erfc(v^|Im( 7 )|). 
Consequently, the investigation of the eigenvalue statistics of GinOE reduces to the study 

f y[Gin]_ 

4.2 The Range of Mahler measure 

The Mahler measure of a polynomial 

N N 

f(x) = J2<inX N - n = a l[(x- 7n ) (4.1) 

71—0 71 — 1 

is given by 

H(J) = exp ( f log \f(e 2 ^ e )\ dd) = \ao\ TT max {^ W}- 

U° > n=l 

The second equality comes from Jensen's formula. Mahler measure arises in a number of 
contexts, i.e. ergodic theory, potential theory and Diophantine geometry and approximation 
(a good reference covering the many aspects of Mahler measure is [TO])- One problem in 
the context of the geometry of numbers is to estimate the number of degree N integer 
polynomials with Mahler measure bounded by T > as T — * 00. Chern and Vaaler 
produced such an estimate in [3] using the general principal that the number of lattice 
points in a 'reasonable' domain in R JV+1 is roughly equal to the volume (Lebesgue measure) 
of the domain. That is, if we identify degree N polynomials with their vector of coefficients 
in R JV+1 and use the approximation 

#{g(x) G Z[x] : degg = N,fi{g) < T} « X N+1 {g(x) G R[x] : degg = N,u.{g) < T} 

= T N+1 X N+1 {g(x) e R[x] : deg 5 = N,n{g) < 1}, 

then the main term in the asymptotic estimate Chern and Vaaler were interested in can be 
expressed in terms of the volume of the degree N star body of Mahler measure, 



Un = {g{x) e R[x] : degg = N^{g) < 1}. 



(4.2) 
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The volume of this set is given by 

— Z w (4.3) 
N+ 1 v ; 

where "0(7) = max{l, |7|} _A,_1 8J- Consequently, the volume of the star body which leads 
to an asymptotic estimate of interest in Diophantine geometry is essentially equals the 
partition function for the random point process . 



4.3 The Range of Other Multiplicative Functions on Polynomials 

We can generalize Mahler measure by replacing the function 7 1— > max{l, |"y|} with other 
functions of 7. Given a continuous function cj> : C — > (0, 00) which satisfies the asymptotic 
formula, 

0(7) ~ H as | 7 | -> 00, (4.4) 
we define the function $ : C[x] — ► [0, 00) by 

AT iV 

* : ao JJ (a; - 7„) h-> |oo| H 0(7") ■ 

n—1 n—1 

The function $ is known as a multiplicative distance function (so named because it is a 
distance functions in the sense of the geometry of numbers on finite dimensional subspaces 
of C[x]). The asymptotic condition in (|4.4| ensures that <E> is continuous on finite dimensional 
subspaces of C[x] (one of the axioms of a distance functions). 
We define the degree N starbody of $ in analogy with (|4.2[) : 

U N {<S>) = {g(x) e R[x] : degg = AT, $( 3 ) < 1}. 

In this situation the volume of Ujv^) equals 

A Ar+1 (U A r($)) = ^- T zM, where </>( 7 ) = fa)-"- 1 . 

We may discover Further information about the range of values $ takes on degree N 
real polynomials by considering the point process on Xjq corresponding to the root function 
"0(7) = 4>{l)~ a f° r cr > N. The partition function of J™ is therefore a function of a, 
and ZM(o-) is known as the degree N moment function of $. In fact, we may extend the 
domain of a real moment functions to a function of a complex variable s on the half-plane 
Re(s) > N. Moreover, in this domain Z^(s) is analytic. Any analytic continuation beyond 
this half-plane gives information about the range of values of $ which may not be realizable 
by other methods. For instance, when ^(7) = max{l, |7|} _s , Chern and Vaaler discovered 
that (s) has an analytic continuation to a rational function of s with rational coefficients 
and poles at positive integers < A^ [5] . Similar results have been found for moment functions 
of other multiplicative distance functions; see [22] . 



5 Correlation Measures and Functions 

Suppose I and m are non-negative integers, not both equal to 0. Then, given a function 
/ : x Cf -> C we define the function F f : X N — > C by 

*xo= E /( x < z )> 

{x,z}C{ 
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where the sum is over all (x, z) e l' x C™ such that {x, z} C £. We take an empty sum to 
equal 0, and thus if f G X LM with L < I or M < m, then F f (}) = 0. 

Given a point process CP on Xjy, if there exists a measure /9^ m on M. e x CJf* such that for 
all Borel measurable functions /, 

f /(x,z)dp*, m (x,z)= / fXOdTU), (5-1) 

then we call pi >m the m)- correlation measure of CP. Furthermore, if pi >m has a density 
with respect to Xi x A2 m then we will call this density the £, m— correlation function of CP 
and denote it by Ri >m . By convention we will take i?o,o to be the constant 1. 

Proposition 1. — — 2lfH — -(x, z) equals 
d(vt x v 2m ) 



2 m 
Z" 7 



^ (L-i)l(M-m)l J J \ A ( x V a , z VP)\ d VL-£{a)dv 2 (M-m){P), 



(L,M) 

L>e.M> 



where x V ex £ M. L is the vector formed by concatenating the vectors xel and a £ M. L ' 
(and similarly for zV/3 £ C M ). 

Corollary 2. If v = wX, then Ri m (x, z) equals 

(L - iy.^M 1 - m)\2 M - m J J n L,M(^o i} zVf3)dX L ^(a)dX 2(M - m )((3). 



(L,M) 
L>£,M>r 



Jl-1 rM-n 



Proof of Proposition^ Assume that L > I and M > m, from (|3.1|) and (|3.2j) . 

/ m)dnA0 = JTM\ i I Ff{<x,P}\H<x,0)\dM<*)dv2M(P), (5.2) 
Jx L , M L\M\J rL J c m 

The function (a,/3) i— ► |A(a,/3)| is invariant under any permutation of the coordinates 
of a. and /3, since such a permutation merely permutes the columns of the Vandermonde 
matrix. Similarly, replacing any of the coordinates of (3 with their complex conjugates 
merely transposes pairs of columns of the Vandermonde matrix. That is, if (cc,/3) and 
(a',P') are elements in R L x Cf such that {a, [3} = {a',/3 1 }, then |A(e*,/3)| = \A(a',/3% 
Moreover, if any of the are real, then the Vandermonde matrix has two identical columns 
and is therefore zero. We may thus replace the domain of integration on the right hand 
side of (|5.2| with R L x C M . In fact, we may assume that the domain of integration on 
the right hand side is over the subset of M. L x C M consisting of those vectors with distinct 
coordinates. 
Now, 

F f {cx,/3}= J2 ^ x ' z )' 

{x,z}C{a,/3} 

Assuming that the coordinates of (ot,/3) are distinct, and if (x, z) G Mr x C™ is such that 
{x, z} C {a, /3} we may find a vector (a, b) £ R L ~ £ x C M ~ m such that (x V a, z V b) is given 
by permuting the coordinates of (a,/3). Clearly |A(a,/3)| = | A(x V a, z V b)| and 

dv L (a) dv%M(0) = dvt(x) dv L -t(&) dv m (z) diy M -m{b). 
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These observations together with an application of Fubini's Theorem imply that (|5.2[) can 
be written as 



Now, it is easily seen that there are 



|A(xVa,zVb)| dis L -i(a) dv M - m (b) > dv e (x) dv m (z). 



2-( L U M )ml 



in 



vectors corresponding to each {x, z}, and thus we find 



{ T 



|A(xVa,zVb)| dv L -t(a) du M - m (b) > dfe(x) dv m {z). 



{L-t)\[M-m)\ 
The proposition follows since 

JX (L,M) L - M 

From here forward, and unless otherwise stated, (I, m) will represent an ordered pair of 
non-negative integers such that £ + 2m < N. 

6 A Matrix Kernel for Point Processes on Xn 

From here forward we will assume that N is even. 

Let J' 1 ' be the point process on Xn associated to the weight w, and as before let v = 
v\ +1^2- We define the operator e„ on the Hilbert space L 2 {y) by 



9(y) sgn(y - 7) dvi (y) if 7 e R, 



ig(l) sgn(Im(7)) if 7 € C*, 

and we use this to define the skew-symmetric bilinear form on L 2 (v) given by 

(g\h)u ■= / g{l)zvh(l) ~ Zv9{i)h{i) du(<y). 



If v = wX then (g\h) v — {g\h)\ where, for instance, ^1(7) := h("f)w("j). 

Theorem 3. Let q = {90(7)1 92(7), ■ ■ ■ , 9jv-i(t)} Q R[t] be such that each q n is monic and 
deg q n = n. Then, 

Z v = PfU£, (6.1) 

the Pfaffian of Uq, where 

U^:=[{q n \q n ,) v ]; n, n' = 0, 1, . . . , N - 1. 
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We will call q a complete family of monic polynomials. 

Remark. It is at this point that it is necessary that N be even, since the Pfaffian is only 
defined for antisymmetric matrices with an even number of rows and columns. A similar 
formula to (|6.ip exists for Z v in the case when N is odd; see [21]. However, we have not 
pursued the subsequent analysis necessary to recover the correlation functions in this case. 

Theorem[3] follows from results proved in [2T]. In fact, [5T] gives a formula for the average 
of a multiplicative class function over the point process on Xn determined by the weight 
function of the real Ginibre ensemble. However, the combinatorics necessary to arrive at 
such averages is independent of any specific feature of Gin(7) and the measure on C specified 
by Gin(7)dA(7) can formally be replaced by any measure v. 

In order to express the correlation functions for the point process y w associated to the 
weight w we will define rj to be a measure on C given by a linear combination of point 
masses, and then use the definition of the partition function and properties of Pfaffians to 
expand both sides of the equation 

= Pf lH A +") (6.2) 

The coefficients in the linear combination defining r\ appear again in terms on both sides of 
the expanded equation, and after identifying like coefficients on both sides of the expanded 
equation we will be able to read off a closed form for the correlation functions in terms of 
the Pfaffian of a matrix whose entries depend on a 2 x 2 matrix kernel. 

In order to define the 2x2 matrix kernel for f w , we let q be a complete family of monic 
polynomials, and we define 



«n(7)==9n(7M7) n = 0,l,...,JV-l. (6.3) 



We then define 



N-l N-l 



Syfr.Y) :=2^ J2 

(7) ^Qn'(l), 

n=0 n'=0 

where we define fi n ,n' to be the n, n' entry of (U™ )~ T . Similarly we define, 

N-l N-l 

(7) e AQV (7') 

n=0 n'=0 

and 

N-l N-l 

DSn^,! 1 ) : = 2 ^2 ^2 AW9n(7)$n'(7')- 
n=0 n'=0 

Remark. The functions Sn,ISn and DSn can be shown to be independent of the family 
q. By setting q to be skew-orthogonal with respect to the bilinear form (-|-)a we arrive at 
particularly simple representations for these expressions. 

Finally, in order to define the matrix kernel CP" 1 we define the function £ : C 2 — > {— h, 0, ^} 
by f 

-sgn(7-V) if7,7'eM; 



£(7,7') := 



otherwise. 
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The matrix kernel of CP™ is then given by 



DS N ( in ') S N ( in ') 
-S N (y n ) IS N ( % ~/') + £(7, V) 



(6.4) 



Remark. The explicit A^-dependence of .Kjv and its constituents is traditional, since one is 
often interested in the TV — > 00 asymptotics of Kjy. 



7 Correlation Functions in Terms of the Matrix Kernel 

We may state one of the main result of this manuscript. 
Theorem 4. The l,m -correlation function of? w is given by 

i^, m (x,z)=Pf 



K N {xj,Xj>) K N {xj,z k f) 
K N (z k ,Xj>) K N (z k ,z k >) 



j,j' = l,2,...,£, 
k,k' = 1,2 



Remark. The matrix on the right hand side of this expression is composed of 2 x 2 blocks, 
so that, for instance, the first row of 2 x 2 blocks is given by 

[K N (xi 1 xi) ... K N (xi,xt) K N {xi,z 1 ) ... K N (xi,z m )]. 

We define the measure 5 on C to be the measure with unit point mass at 0. Given U 
real numbers Xi, X2, ■ ■ ■ , xu and V non-real complex numbers, zi, Z2, . ■ ■ , Zv, we define the 
measure r\ by 

u v 
df l(l) : = X! Qu ~ Xu ) + ( d6 ^ ~ Zv ^> + d6(j -z^)), 

u—l v—1 

where 01, a%, . . . , au and 61, hi, ■ ■ ■ , by are indeterminants. It does no harm to assume that 
U and V are both greater than TV. By reordering and renaming the x, z, a and b we will 
also write 

T 

d vh) -=^2^(^(1 -yt), 
t=i 

where we define 



dd~h~y) if ye 

d8{n — y) + dS(j — y) if y 6 



Clearly T = U + V. 

As we alluded to previously, the proof of Theorem 0] relies on expanding Z w ^ x+V ' ) /Z w 
in two different ways and then equating the coefficients of certain products of c%, 02, . . . , cy. 
One of the expansions comes from Theorem|31 while the other comes directly 

from the definition of the partition function. 



Proposition 5. 



Pf(j+ [V^K N (yt,yf)]); t,t' = 1,2,...,T, (7.1) 
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where 3 is defined to be the 2T x 2T matrix consisting of 2 x 2 blocks given by 

tt=l,2,...,T. 





' 1" 




6 t ,t> 


-1 





Remark. The Pfaffian which appears in (|7.ip is an example of a Fredholm Pfaffian. This is 
the Pfaffian formulation of the notion of a Fredholm determinant and is discussed in l 20] . 

We defer the proof of Proposition [5] and Proposition [6] until Section [9] 

For each (£, m) and (L, M) we define the £, m, L, M-partial correlation function of 3™ to 

be Rt>, m ,L,M :l'xC m ^ [0,oo) where 



l,m,L,M (x, z) : = 
1 



{L-£)\{M -m)\2 M - r 



/ / fi iiM (xV q,zV /3)dX L -e(a)d\ 2{ 

M-m) 63). 



When I = m = 0, we take x V a. = a and z V (3 = (3, so that Rvi.o,l,m is a constant equal 
to 5 ) l,m(^l,a-/)- The partial correlation functions are related to the correlations functions 
by the formula 

i^, m (x,z) = — ^2 ^,m,i,M(x,z), (7.2) 

£>f ,M>m 

and thus the partial correlation functions are one path to the correlation functions. Equa- 
tion (|7.2p is still valid when I = m = 0, though the constituent correlation 'functions' are 
actually constants. 

The partial correlation functions of CPl Gm ] of the special forms Rl,m,o.m and Ro,N,o.n 
have been studied by Akemann and Kanzieper in [1] and [17j . For general weight function 
w, the partial correlation functions of the form RN,o.n,a are (up to normalization) equal to 
the correlation functions of the — 1 Hermitian ensemble with weight w. This connection 
will be exploited in [XJ 

Before stating the next proposition we need a bit of notation. Given non-negative in- 
tegers n and W, we define 3^ to be the set of increasing functions from {1, 2, . . . , n} into 
{1,2,..., W}. Clearly if W < n then 3^ is empty. Given a vector a £ C w and an element 
u £ 3jf , we define the vector a u £ C™ by a u = {a u (i), a u(2 ), . . . , a u( „)}. 

Proposition 6. For each pair (L, M), 

1 



L\M\2 M 



n LM (a,p) d(X + r)) L (a) d(X + n) 2M {P) = 

L M f I m *, 

EE E E 1 IT «u(j) II & «(fc) [^m,L,M(x u ,Z B ) 
fcO rn=0ue3^ D£3^ ^ j=l fe=l J 

Remark. We will use the convention that 

o o 

E II a «W = J2 n & »w = L 

ue3, f ,'i=i veoX fc=1 



This will allow us to keep from having to deal with the pathological correlation 'functions' 
Rq,o,l,m and i?o,o separately. 
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Proposition 7. Suppose K : {1, 2, . . . , T} X {1, 2, . . . , T} 

K(t,t')=-K(t',t) T , 



is such that 



and define K to be the 2T x 2T block antisymmetric matrix whose t, tf entry is given by 
K{t,t'). Then, 



Pf[J + K] = l + E J2 PfK '' 



(7.3) 



where for each t G K t is £/ie 25 x 25* antisymmetric matrix given by 
K i =[^(t(s),t( S '))]; a.a 7 = 1,2,..., 5. 

Proof. This is a special case of the formula for the Pfafhan of the sum of two antisymmetric 
matrices. See [22] or [23] for a proof. □ 

Using these lemmas we may complete the proof of Theorem 0J First notice that Propo- 
sition [6] and (|T.2[) imply that 



, I rn . 

E E E \ II a «0') Ii b m f#«,m(Xu,Zt>) 
(2,m) ue3? oea^ ^ i =1 fc=i J 



(7.4) 



This follows by summing both sides of the expression in Proposition [7] over all (L, M) and 
then reorganizing the sums over (L,M), £ and m. 

Given u £ 3f and <E 3^ we define (u: d) to be the element in given by 



(uV t>)(n) 



u(n) if n < ^ 

^ + o(n) if n > 1 



Notice that each t G 3^ is equal to (u V d) for some u£ 3^ and € 3^ with £ + rn = S. 
(This does not preclude the possibility that either U or V equals 0, in which case t = D or 
t = u) . It follows that we can rewrite (17.311 as 



pf(j + K)= e EE pfK (UVD) 

l + m<T 

If we set K t ,t'(t,t') = ^c t c t > K N (y tl y t '), then 



Pf K 



(uvo) 



n a »o) n 6 »w r pf 

3=1 fc=l 



^Ar(^u(fc))^u(fc')) K N (x u (k), Zo(n')) 

K N { x u(k')) K N ( 



where k and k' are indices that run from 1 to £ and n and n' are indices which run from 1 
to m. Thus, Pf(J + K) equals 



E E E n«»o)IW) pf 



-^Ar(z D (n): ^u(fc')) ^( Z n(n)' z o(n')). 



, (7.5) 



and Theorem 2] follows by equating the coefficients of Ylj=i a j IlfcLi ^fc m (|7.4p and (|7.5[) . 
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8 Limiting Correlation Functions for the Real Ginibre 
Ensemble 

We now turn to the large N asymptotics of the matrix kernel for the Ginibre ensemble of 
real matrices. In fact, we maintain our restriction to the case where N = 2M is even, and 
consider the asymptotics of Kim as M — > oo. Throughout this section we will take cj> to be 
the function given by 



0( 7 ) = exp(- 7 2 /4 - 7 2 /4) v /erfc(\/2|Im(7)|). 



Notice that, since erfc(O) = 1, when 7 € E this reduces to exp(— 7 2 /2). 

The skew-orthogonal polynomials for this weight are reported in [T^] to be 



7r 2m (7) = 7 2m ^+1(7) = 7 2m+1 - 2m 7 2 "'- 1 



with normalization ("K2m\'K'2m+i) u = 2V27r(2m)!. A detailed account of the derivation of 
these skew-orthogonal polynomials will appear in |llj . 

The skew-orthogonality of these polynomials and formulas from Section [6] imply that 



Q f '\ 1 7r 2m(7)eA7T2m+l(7') ~ 7r 2m+l(7)eA7T2m(7') / Q , s 

W7 ' 7 } = ^ ' (8 } 



m=0 
M-l 



nC / M 1 \- 7r2m(7)71-2m+l(7)-7r2m+l(7)7T2m(7) , Q x 

Db 2 M(J,J ) = —f= > 77 — T , (8.2) 



m=0 
M-l 



ro / A 1 eA7T2m(7)eA7r2m+l(7 ) - eA7T2m+l(7)eA7r2m(7 ) /o o \ 

^2M(7>7 ) = "7= >, 7^" ■ ( 8 -^) 

V27T ^— ' (2m)! 

v m— v 7 

The correlation functions are all of the form Pf K for an appropriate matrix K whose 
entries are given in terms of (|8.1[) . (|8.2p and (|8.3[) . If D is a square matrix such that the 
product DKD T makes sense, then Pf (DKD T ) = Pf K • detD. And thus, if detD = 1 we 
have Pf K = Pf(DKD T ). That is, we may alter (and potentially simplify) the presentation 
of the Pfafhan representation of the correlation functions by modifying K in this manner 
by a matrix with determinant 1. When D is diagonal, the process of modifying K by D 
preserves the block structure of K. That is, the effect of modifying K by D affects changes 
at the kernel level and the correlation functions can be represented as the Pfaffian of a block 
matrix (c/. Theorem 2]) with respect to a new matrix kernel dependent on Kjy and 
D. This will allow us to write the correlation functions of the real Ginibre ensemble in the 
simplest manner possible by 'factoring' unnecessary terms out of the kernel. 

It will be convenient to define cm, sm and eM to be the degree 2M—2 Taylor polynomials 
for cosh, sinh and exp respectively. Explicitly, 

rfrn. M-l 2m-l 
c M{l) ■■= V *Af(7) : = V 77, 7TT 

(2m)! '-^ (2m— 1)! 

m=0 v ' m=l v ' 

and 

2M-2 „, 
7- 

ml 

m=0 
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Theorem 8. The I, m-correlation function of the real Ginibre ensemble of 2M x 2M ma- 
trices is given by 



i^, m (x,z) = Pf 



K 2M {Xj,Xj>) K 2M (Xj,Z k >) 
K2M(zk,Xj' ) K-2M (Zk , Zk> ) 



j,f ' = 1,2,..., £; 
k , k — 1,2,..., Tn^ 



where 

y-S 2M {l',l) IS 2M {l,l') + £{l,l') 
is given as follows. Let x, x' e R and z, z' e C, i/ien: 

i. T/ie entries in the real/real kernel, K 2 m{x,x'), are given by 



S 2M {x,x') = 

= -z 2 /2 



e 2 



2tt 

>: v v ^ Im ^)(2.U-2,! 



e 11 eu(xx') + tm{x,x'), where 



e~ z jA I p 

r M (z,x) = — r^v/erfc (v2Im(z)) 



2M-3/2 



and 7 is i/ie (oraer incomplete gamma function; 
l\ c / / 



• DS 2M (x,x ) 

• IS 2M {x,x') 



2tt 



-{x' - x)e xx e M {xx'); 



7 f 
— sgn(x') / 

^ Jo 



x' 2 /2 e -t 



, - , c M {xV2t)dt - _ , . , 

2. The entries in the complex/ complex kernel, K 2 m{z, z'), are given by 

• S 2M {z,z') = lC 2{Z * ) (z' - z) yl / crfc (V2Im(z)) crfc (V2lm(z'))e- Zz ' e M {zz'); 
V 2?r 



-x' 2 /2 r 

r^sgn(x) / 
Jo 



* 2 / 2 e^ 



c M (x'V2t)dt. 



DS 2M {z,z') = -;=—{*! ~ z)^eTfc(Vzlm(z))eTfc(V2Im(z>))e- zz 'e M {zz'y, 

z) y'crfc (V2Im(z))erfc (v^Im(^)) e "^' e m(zz')- 



e -|(z-z') 2 
• IS 2 m(z,z / ) = — — (2 



5. TTie entries in the real/complex kernel, K 2 m{x,z), are given by 

zc ^ ^ / 

• S 2 m(x,z) = = — (2 — z)y erfc (v^Im(z))e~ X2: eM(£z); 



'2tt 

_ g-^(2:-z) 2 ; 

• SW^z) = 1= — verfc (V2Im(z))e~ xz e M (xz) +r M ( 

V27T 

• DS 2M (x,z) 



z,x . 



2ir 



(z - x)Jer{c (\/2Im(z))e xz e M (xz); 



IS 2 m(x,z) = j= — \/erk(V2lm(z))e~ xz eM(xz) -irM{z,x) 



2n 
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Theorem [8] allows us to derive the M 



limit of K- 



2M- 



Corollary 9 (Limit at the origin). Let x and x' be real numbers, and suppose z and z' are 
complex numbers in the open upper half plane. We define K = lim Kim- Then, the limit 

M—too 

exists, and 



1. The limiting real/real kernel, K(x,x'), is given by 



K(x,x') 



-h={x' — x)e 



-±(x-x'f 



=e 2 



\{x-x'f 



1 e -|(z-z') 2 

V27T 



| sgn(x - x') erfc 



2. The limiting complex/complex kernel, K(z,z'), is given by 
K(z,z') = -^Jeiic (V2lm(z)) erfc (V2lm(z')) 



(z' - z)e~^ z ~ z "> 2 i(z' - z)e~^ z - z "> 2 
i(z' - z)e-5(--^') 2 _ ^ e -±(---') ; 



3. The limiting real/complex kernel, K(x,z), is given by 



K(x, z) 



'2tt 



erfc (v / 2Im(z)) 



— t\p. 2 ( X 



(z — x)e 



(x-*) a 



i(z — x)e z( x z ^> 
—le ? y ?■> 



Remark. Observe that all blocks of the kernel are invariant with respect to real shifts. That 
is, if c G K and 7, 7' are in C then 

K(-f + c, 7' + c) = if (7, 7'). 



8.1 In the Bulk 

The circular law for N x N matrices with i.i.d Gaussian entries says that, when normalized 
by N^ 1 / 2 the density of eigenvalues becomes uniform on the unit disk as N — + 00 (See [15] 
for a proof of this fact when the entries are i.i.d. Gaussian, and [2] and [24] for more general 
results). This gives us the appropriate scaling when considering the matrix kernel in the 
bulk. Specifically, in this section we will be interested in the large M limit of K2m(uV'2M + 
s, u\JlM + s') where u is a point in the open unit disk, and s and s' are complex numbers. 

When u is real we expect that the limiting kernel under this scaling should yield K(s, s'); 
indeed this is the case. When u is nonreal a different kernel arises. 

Theorem 10. Let —1 < u < 1 be a real number, let r%, rj, . . . , r< G 1 and s\, S2, ■ ■ ■ , s m be 
in the open upper half plane. Set, 



Xj = u\2M + rj j = 1,2,..., £ 
Then, 



Zf~ = uV2M + Sfc k — 1,2,..., to. 



M 



lim i?£ iJn (x,z) = Pf 



K {rj,r f ) K( rj ,s k >) 
K(sk,rji) K(sk,s k ') 



j,j' = 1,2,... J; 
k, k' = 1, 2, . . . , to, 



where K is given as in Corollary^ 
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Theorem 11. Let u be in the open upper half plane such that \u\ < 1 and suppose 
si, s 2 , . . . , s m G C. Set, 



Zk = tiv 2M + Sk 



Then. 



lim Ro tm (— ,z) = det 



exp 



1,2,. 



\Sk' 



SkSk' 



k,k' = l 



Remark. The limiting correlation functions in the complex bulk are invariant with respect 
to any complex shift. 

Remark. The function 



i / |s| 2 |s' 12 



M~-exp(-— - — + .s 

is, up to a factor of 1/2, the limiting (scalar) kernel of the complex Ginibre ensemble. Thus, 
the limiting correlation functions in the bulk of the real Ginibre ensemble off the real line 
is almost identical to the limiting correlation functions in the bulk of the complex Ginibre 
ensemble. See Ginibre's original paper [14], or [T9J, Section 15.1], for the derivation of the 
finite N correlation functions for the complex Ginibre complex ensemble. We derive the large 
N asymptotics of the correlation functions of Ginibre's complex ensemble in Appendix [Cl 



8.2 At the Edge 

At the edge of the spectrum new kernels emerge. 

Theorem 12. Let u = ±1, let r\, r 2 , ■ ■ ■ , Tt € M and si, s 2 , ■ ■ ■ , s m be in the open upper half 
plane. Set 



x-j =uV2M + rj j = 1,2,. 
Then, 



■ A 



and 



z k = uV2M + s k 



k = 1,2,. 



M 



lim i?£ jm (x, z) = Pf 



K cdgc {r ,r r ) K eds e(rj,s k >) 
Kei ge {s k ,r r ) K cdgc (s k :Sk') 



J,/ -1,2,. 
k,k' = 1,2, 



■A 



whe 



DS edgc (7, 7') Sedge (7, 7') 
-Scdgch',7) ISedgeil,! 1 ) 



and: 



1. The real/real kernel at the real edge, K cdgc (r,r'), is given by 

• Sedgc(r-,r') = -l= e -^ r ~ r "> 2 erfc (u ^p- ^j + erfc(-ur'); 



• DS t 



IS edgc (r, r') = - sgn(r - r') erfc 



V2 
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2. The complex/ complex kernel at the real edge, K e d ge (s,s'), is given by 

• <Sedge( S 7 S ) 



2V27T 



rfc (\/2Im(s)) crfc (V2Im(s')) 



x(s' - s ) e -^ s - s "> 2 erfc I u 



(* + *') 
V2 



D5edge(s,s') = -^=\/erfc (\/2Im(s)) erfc (\/2Im(s')) 
2\/27r v 



IS c dgc{s, S ) 



1 



:( S '- S )e-^ S ') 2 erfc(u^±|^ 



2V2tt 



erfc (V2Im(s)) crfc (\/2Im(s')) 

x(s' - s)e _ ^ s_s ') 2 erfc 



+ s' 



5. T/ie real/complex kernel at the real edge, K e ^ ge (r, s), is given by 
• S cdgc (r,s) = —!— e -Ur-sf^ CTic (V2Im(s))(s - r) crfc 



Sedge ( s 7 7") 



2V27T 
1 



(r + s) 
' V2 



-|(r- S )' 



2V27T 



e 2 



erfc (-\/2ini(s)) erfc 



(r + s) 
1 V2 



+- — =e s erfc(— ur); 



£>5 edge (r, s) = -l=e-^ r - a ^ ^crfc (v^Im(s)) (s - r) crfc 



(r + s) 



• /Sedge (r,s) 



-t _ l 



(r-s) 5 



2V27T 



erfc (v / 2Im(s)) erfc I u 



(r + s) 
V V2 



e s erfc(— ur). 



4^ 



Remark. The kernel when u = — 1 is the image of the kernel at u — 1 under the involution 
on the closed upper half plane given byz^ —z. 
At the complex edge we have the following: 

Theorem 13. Let u be in the open upper half plane such that \u\ — 1 and suppose 
Si,S2, ■ ■ ■ , s m e C. Set, 

z k = W2M + s k k = 1, 2, . . . , m. 

Then, 

— I - XT'" 
SfeU + Sfe/ti 



M 



lim i? ,m(^7 z ) = det 



1 

— exp 



|s fc | 2 \s k '\ 2 



+ s/cSfc' erfc 



V2 



k,k' = l 



Remark. The kernel at the complex edge are identical to that of the kernel at the edge of 
the complex Ginibre complex ensemble. 
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9 Proofs 

9.1 The Proofs of Proposition [5] and Proposition IA.31 

In the case of the real asymmetric ensembles Y = C and 6=1. In the case of the Hcrmitian 
ensembles Y = M and b = y/]3. We start with 

{qn\q n ') w (\ + r,) = {q n \q n ')x+n = (q n (i) ^x+rfiwii) - ex+^qnii) q n >{i)) d(X + 77) (-y) ■ 

An easy calculation reveals that this is equal to 



\Qn I qn' ) w X + T^2ct(q n (yt)exqn'(yt) - qn'(yt)exq n {y u )) 



T T 



t X! X! c t c t'qn{yt)q n >{yt>) £{yt,vt')- 



t=i t'=i 

Next we define A to be the Nb x 2T matrix given by 



A := 



q n {yt 



b V b 

and the 2T x 2T matrix B. 



txq n (yt) 



n = 0,l,...,Nb-l and f=l,2,...T, 



B 



J + 



y/ctcr £(yt,yt>) 




t,t' = 1,2,..., T. 



We define C to be the Nb x Nh matrix given by C = (U^ A ) T ; the n, n' entry of C is \i T 
A bit of matrix algebra reveals that 

Z w(x+ n ) ^ Pf ( C -T - ABA T ) _ Pf(B- T - A T BA) 



Pf(C- T ) 



Pf(B- T ) 



(9.1) 



where the second equality comes from the Pfaffian Cauchy-Binet formula (see [Bj . Notice 
that 



B 



tj = 1,2,..., T 



-J- ° 

[0 y/ctct' £(yt,w)_ 

And from Proposition [71 Pf(B) T = Pf(— J) = (— 1) T . A bit more matrix algebra reveals, 



A'CA = 



y/ctCf DS N (y t ,yt>) y/c t ct> S N {y t ,y t >) 
-y/ctct' S N (y t ',y t ) y/ctCf IS N (y u y t >)_ 



t,t' = 1,2,... T. 



Using these facts and simplifying (|9.ip we find 

yw(X+r]) 

zw = (-i) T Pf ( - J - WwKN{vt,w)]) ; t, t' = 1, 2, . . . , t, 

and the Lemma follows by using the fact that if E is an antisymmetric 2T x 2T matrix, 
then Pf(-E) = (-l) T Pf(E). 
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9.2 The Proofs of Proposition [6] and Proposition IA.41 

We start with 

' ' ' n L M(a,P)d(\ + i 1 ) L (a)d(\ + r ] ) 2M ((3). (9.2) 



L\M\2 M 



Notice that in the case of the Hermitian ensembles that this is equal to Z w ^ x+r) ^ when L = N 
and M = 0, and the proof of Proposition DOl follows from the proof recorded here by setting 
every instance of M to 0. 
First we write 

L L 

d{\ + v)t( a ) = IT (dXi(aj) + dr)i(aj)) = dr}i(a u )dX L -e{a u >), (9.3) 

3=1 i=0 u£3f 

where given t £ 3^ , we define t' to be the unique element in 3^_ n whose range is disjoint 
from t. Notice that since u' appears in the summand on the right hand side of (|9.3|) . the 
inner sum is not actually empty when £ = 0; in this situation the summand is equal to 
dAi(a). 

Similarly, 

M 

(M-m) (PA 

m=0 o£3 M 



Thus, JOJ equals 

L M Iff 

EE E E L\M\2 M I I ^LM( a >P) d Ve(atu)d\ L -e(a u >)dri2 m (f3 )d\2(M- m )(P >)- 

We can relabel the a and [3 in the integrand in any manner we wish, and in particular we 
may make the integrand independent of u and 0. In particular, if we set i £ 3^ to be the 
identity function on {1, 2, . . . , n}, and since the cardinality of 3^ is CY) , we find that (|9.2[) 
is equal to 

l M 



EE 



^ l\(L - t)\m\(M - m)\2 M 
x J f {^f e J ^LM( a ,0)dve{ai)dij2 m {Pi)^dX L -e(ot i ,)dX 2 (M^ m )((3 i ,). (9.4) 

H r — < C M - m 

Now, 

t IV 
d-q t (oLi) = _Q??i(ai) = If }^a u dS(aj - x u ). 

We may exchange the sum and the integral on the right hand side of this expressions by 
using the set, ^ °f au functions from {1, 2, ... ,£} into {1, 2, ... , U}. Specifically, 



dm(<*i) = E \ II a "C0 dS ( a J ~ Mi)) f> 

..r-f*U y a— 1 ' 
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and similarly, 

l, A' fe=l 



Thus, 



f2z,,M(«, /3) drjiiati) dr?2m(A) 
ups^oes, 1 : M=i fc=i J 



(9.5) 



Notice that if u or is not one-to-one then |A(x u V on, z D V /3 t ) j = 0. We may consequently 
replace the sums over §Y and $\ with their respective subsets of one-to-one functions. 
Moreover, since &l,m is symmetric in the coordinates of each of its arguments, we may 
replace each one-to-one function in these sums with the increasing function with the same 
range so long as we compensate by multiplying by l\ and ml. Proposition [6] follows from 
the definition of Ri, m ,L,M by substituting (|9.5[) into (|9.4[) . Proposition follows from the 
fact that R n = R n ,o,N,o/Z w . 

9.3 The Proofs of Theorem M and Corollary M 

It shall be convenient to introduce the following variants of S2M and DS2M'- 

£ , ,x 7T2m(7)eA7f2m+l(7') - ^2m+l{l)^2m{i) 

2m)! 

m=o v ; 

and 

775 / is 7r 2m (7)7r 2m+ i(7') - 7r 2m +i(7)7r 2m (7') 

^-^ (2m ! 

The following lemma gives a closed form for these functions. 

Lemma 9.1. Let x be a real number, and suppose z and z' are complex numbers. 

1. S 2M {z,x) = <p(x)e M (zx) + (2M-2)! 8gn ^ z2tf ' 1 -jlM- -, — 

2. DS 2M (z,z / ) = {z> - z)e M {zz'). 

Proof. First we compute eyK2m and eyK2m+i- We start by noticing, 

1 f°° 1 f x I f°° 

£\9{x)=2 9(y)sgn(y-x)dy = -- g(y)dy+- g(y)dy. 
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When g(y) = e v ^ 2 y n , we may evaluate the latter two integrals in terms of the incomplete 
gamma functions. 



2(n-l)/2p^n+l x^_\ 

We immediately conclude that 



2 (n-i)/a ggn^) . 7 (™±I, ofj if n is even; 

if n is odd. 



e\-K2 m {x) = -2 ro ^ 2 sgn(a;) • 7 fm + ^, , 



and 



CA^m+i^) = 2'' 



r m + 1, 



m,r r 



(9.6) 
(9.7) 



where in the second equality we used the fact that T(a + l,x) = aT(a, x) + x a e 
Using (|9.6p and (|9.7p . we may write 



S 2 m{z,x) = cj>(x)c M {zx) 

E 72^r z2m+l7 i "' 

m=0 ^ 



„2\ M ~ 1 2 m ^ 1 / 2 



1 a; 

2 '"2 / ^ (2m- 1)' 

7 m=l v ' 



E 



z 2m X 7[ TO 



1 x 2 

2' y 



Next, we use the fact that 

j(a + 1, x) = 07(0, x) — x a e~ 
so that the second sum in this expression becomes 



(9.8) 



M ~ 1 o(m-l)-l/2 



E 



r 2(m-l) + l r 



^ (2(m-l))! 
Consequently, 



1 a; 



7((m-l) + -,y ) +<(>(x) E 



W_1 z 2m-l| x |2m-l 



Sz M (z,x) = <P(x)(c M (zx) +sgn(x)s M (z\x\)) + 



0M-3/2 
(2M- 2)! 



(i)eM(zi) + 



9JW-3/2 
(2M- 2)! 



sgn^z^-^lM-^, 



(2m- 1)! 



1 x 2 



2' 2 



Turning to -DS* 



A ^ z 2 ™(z' 2m + 1 - 2mz /2m - 1 ) - (z 2m+1 - 2mz 2m - 1 )z' 2m 
DS 2 m(z,z ) = 



E 

m=0 



(2m) 



(2m)! 

M 

E 



m— 1 



(2m- 1)! 



,2m,/2m M-l 2m-l,/2m- 1 



Z) j ^ (2m)! + ^ (2m- I.)! 
1^=0 v ; m=l v ' 

= (z' - z)e M (z'z). 



□ 
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With a closed form for S 2 m and DS 2 m in hand, we are ready to prove Theorem [8] 
Proof of Theorem [3 From Lemma 19.11 we have that 

t>2M\x,x ) = —==b 2M {x,x ) = -== — e M \xx ) + r M (x,x ) 

VZTT s/2lT 



and 

DS 2m {x,x ) = -== — DS 2 m[x,x ) = y== — (a; - x)e M {xx ). 

V27T \J2lT 

Now, 

and therefore, 
and 



x)<f>{x') = e -5(- 2 +-' 2 ) = e -W-*')\-™\ 



g 2 (x X ) f 

S 2M {x,x') = — — e~ xx e M {xx') +r M (x,x'), 

x/2n 



DS 2M (x,x') = = — (x' - x)e~ xx ' e M (xx') 

\/2ir 



The computation of IS 2 m(x,x') is a bit more involved. From (|8 . 3[> . (|9.6[) and (|9.7p . we 

1 A f~} o m / 1 ■r' 2 \ 9 
IS 2M (x,x>) = ^ £ ^ • 7 (m + -, ^) s g n(*'>-e-* / 2 



m=0 

M-l 



2^/tF ^ (2m)! 'V 2' 2 

v m=0 v ' v 

IS 2 m{x,x') is clearly skew-symmetric in its arguments; looking at the first sum in this 
expression we thus find, 

1 M ~ 1 2 rn , , r x ' 2 / 2 



om rx /z 

- T 7^sgn(x')x 2m e' x ' 2 / f^^e^dt 
20F^ o (2m)! J 

= ^e- 2 / 2 sgn(x') / * £ T^T*" 1 ^ ^ 

2v^ 7o V* fr'n ( 2m ) ! 



, m=0 

where on the left hand side we have replaced the lower incomplete gamma function with its 
integral definition, and on the right hand side we have exploited the linearity of the integral. 
The sum on the right hand side of this equation is equal to CM{xV2t), and thus 

IS 2M {x,x') = 

e -x 2 /2 rx' 2 /2 e -t e -x' 2 /2 l-x 2 /2 £ -t 

— -^sgn(V) / —c M (xV2t)dt sgn(a;) / — c M (x W2t) dt. 

W n Jo yt Jo V* 

Turning to the complex/complex entries of K 2 m, if z is assumed to be in the open upper 
half plane then e\ir n (z) = iiT n (z). From this we see that S 2 m(z, z') — i(j){z)(j){z l )DS 2 M{z,z^), 
DS 2M (z,z>) = cP(z)</>(z')DS 2 m(z,z') and IS 2M {z,z') = -<t>{z)<t>{z')DS 2M {z,T?). 
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Next, we define 
Notice that 

t)}{z)'i)}{z )e~ 

and thus 



4>{z)4>{z') = ^(z)^(z')e~^ z ' z "> 2 -y/erfc (y/2 Im(z)) erfc (y/2 lur{z)')e- zz ' . 
Using this and Lemma ISTI! we conclude that 

SWC^O = il>(z)il>(t) — ' (z' - z)^/erfc (V2Im(z)) erfc (%/2Im(z)')e~ z ~'e M (^'), 



/2tt 

i (z -,') 2 

and 



DS 2 m{z, z ') = ip( z )ip(z') e 2 — (z' - z) J eric (V2Im(z)) erfc (^k(z)')e- B 'ejf(z/), 



5 2 Af(z, z') = i>(z)i>(z') C 2( * Z ) {z' - z^crfc (\/2Im(z)) erfc (V2Im(z)') e -"'eAf (zz'). 



/27T 

Lastly we look at the real/complex entries of K 2 m- As in all other cases, 
Db2M\x,z) = DS 2M (x, z), 

V 27T 



and it is easily verified that 

0(x)0(z) = ^(z)e-^ x - z ^^eric(V2lm{z))e- xz . 

Thus, 

g-i(x-z) 2 , — 

DS 2M (x,z) = ip(z) = — (z - x)^Jerk (V2 Im(z))e xz e M (xz). 



/2tt 
Since, 

S 2 m(x,z) = i — L>A 2 m (x, z), S 2 m{z,x) = — ==6 2 m(z, z)- 

V27T V27T 



and 

IS 2M (x 1 z) = -i^=S 2 m(z,x), 

V 27T 

It follows that 



jg— f(x-z)~ / — 

S , 2m(^, z) = ip(z) j= — (z — x) <J eric {v21m(z))e~ xz eM(xz), 



( e -±{x~z) 2 , — _ "I 

S 2M (z,x) = ip(z) < -== — yerfc (V2Im(z))e xz e M (xz) + r M (z,x) \ , 
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and 



I e ~h( x ~ z ) i _ 

IS 2 m(x,z) = -iip(z) < -= — yerfc (\f2lrci(z))e~~ xz eM{xz) + r M (z,x) 



Clearly, ifj(x) = ij)(x') = 1, and thus we find that 



o vCr) 



o 



It follows that, if we define K to be the matrix 



K 



K N (xj,xj>) K N (xj,z k >) 
K N (z k ,Xf) K N (z k ,z k >) 



J,f = 1,2,... J 
k,k' = 1,2 



5 

,m, 



and D to be the diagonal matrix 

D = diag(V'(xi), ^(xT), . . . , ip(xi),tp(xt),ip(zi),ip(zi) . . . ,~ip(z m ),ip(z^)), 

then 

ft, m (x,z) = Pf(DKD). 
But, since ip(~z) = -0(z) _1 , we have that detD = 1, and i?£ >m (x, z) = Pf K as claimed. □ 

Proof of Corollary We first make use of the fact that 

lim eif(z) = e z 

M— >oo 

pointwise on C. This simplifies all terms in the kernel except the IS2M term. It remains to 
show that rufz, x) — ► as M — > 00, and that 



— sgn(x — x 1 ) + lim IS2u(x, x') = — sgn(x — a;') erfc 

2 M— >oo 2 



|ar — a/| 

~7T~ 



(9.9) 



The first of these facts is easily seen by noting that j(M — 1/2, x 2 /2) < T(M — 1/2), and 
by Legendre's duplication formula, 

| z |2Af-i 2 m-3/2 r(M _ 1/2) ^ | z |2M-i 2- M + 3 / 2 V^T(2M - 2) 



T(2M - 1) 



r(2M - i)r(M - 1) 



Thus, 



Vm 



{z,x)\ < e -5^(^ 2 ) J erfc (V2Im(z)) 



I2M-1 



2 M r(A/) : 



and it is easy to see that this goes to as M — > 00, independent of the value of z. 
To establish (19.911 we start with 



I M {x,x') 



—j=CM (xV2t) dt 
yt 



(9.10) 
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Since the terms in cm are all positive, from the Monotone Convergence Theorem, 

• cosh(xv / 2t) dt 

Jo Vt 



I(x,x) := lim Im{%,% 1 ) 

M— j-oo 



v/2 



erf 



\x'\-x 



V2 



(9.11) 



The latter equality follows from the fact that 



o V~t 



sh(aVi) dt 



erf (| + ^/xj - erf - Vi) 



which can be verified via differentiation. 
Now, 

lim IS 2M (x, x') = -1= [c^l 2 sgn(x')I(x, x') - e~ x ' 2 ' 2 sgn(x)I(x', x)\ 

A/-»oo 2y7T L J 

1 



sgn(x') erf 



\x'\+x 



V2 

— sgn(x) erf 



— sgn(a: / ) erf 
x' + \x\ 



V2 



\xf\-x 
+ sgn(a;) erf 



x' - \x\ 
V2 



-— erf 
2 



V2 



It follows that (|9.9[) can be written as 



Isgn^-xO-ier^ — 



= — sgn(x — x') — — sgn(x — x') erf 



\x-x'\ 
V2 



where we have exploited the fact that erf is an odd function. We arrive at the form for 
IS2M stated in the corollary using the fact that erfc = 1 — erf. □ 

9.4 The Proofs of Theorem [10] and Theorem [TTI 

In order to prove Theorems \W\ and [Til it is necessary to investigate the asymptotics of the 
partial sums the exponential function. 

Lemma 9.2. Let u 7^ ±1 be a complex number, and let (u m )m=i be a sequence of complex 
numbers satisfying 

v M =u 2 + (il:T 1/2 ) as M — > 00. 

Then, as M — > 00, 



-2Mv M 



e M (2Mv M ) ~ 1 - 



-2(l-u 2 ) 2M(l-u 2 ) AM 



2ttm 2 (1 - u 2 ) 



In particular, when u is real and < lit] < 1. 



lim e- 2MvM e M (2Mv M ) = 1. 

M->oo 
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Proof. Set v = v M - We start by writing 2Mv = v(2M - 2) + 2v. Thus, 



-2Mv 



e M (2Mv) = cxp -(2M - 2) [v + 



M - 1 



e M (2M - 2) [v + 



M — 1 



We write 



10 = Wm =V + 



M- 1 



Clearly u>m = u 2 + 0(A/ 1 / 2 ). Under this hypothesis, and since w 2 ^ 1 , equations (2.9), 
(2.15) and (1.7) of [5] imply that 



e-( 2M -V w e M ((2M-2)w)~l 



e -2(l-u 2 ) 

2nu 2 (l-u 2 ) 



2M(1- 



AM 



AI 



The second statement of the lemma follows from the fact that if u is real and < \u\ < 1, 
then 



e 2M(l~u 2 ) u 4M _ e 2M(l-u 2 +2IogM) 



and 1 — u 2 + 2 log \u\ is negative when u is in (—1, 1). 
We are ready to prove Theorem [TT1 



□ 



Proof of Theorem \11[ Let u be a point in the open upper half plane with modulus less 
than 1, and suppose s and s' are complex numbers. For all but finitely many values of M, 
z = u\j2M + s and z' — u\j2M + s' are in C*. Thus, in this case, we need only consider 
the asymptotics of the complex/complex kernel under these substitutions. 
We will make use of the fact that if x is a real number, 



erfc(x) 



> (V 

A \2' 



-AM 



Consequently, 



erfc(\/2Im(uV2M + s)) 
Now, by Theorem [51 

DS 2 m(uV2M + s, uV2M + s') 



exp ( - 2Mlm(u) 2 - 2V2Mim(u)Im(s) - Im(s) 2 



^j2lm{u)^Mn 



(9.12) 



'2tt 



'erfc (V2lm(uV2M + s)) erfc (V2lm(uV2M + s')) 
x exp(-2Mw 2 -u\/2M(s + s') - ss'^j e M (2Mu 2 + uV2M(s + s') + ss'^j 
Therefore, by Lemma |9"T2"1 and ([9~T2]) . 



DS 2M (uV2M + s,uV2M + s') ~ ^— *K -^ S - S 'f 

v2tt 



e -4J\fIm(ti) 2 e -2\/2MIm(«)(Im(s)+Im(s')) e -Im(s) 2 -Im(s') 2 

2Im(u)VM7T 

e -2(l-u 2 ) e 2M(l-« 2 ) u 4A/\ 



X 1 



2ttw 2 (1 - u 2 ) 



A I 
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It is easily seen that 



lim • 



-4MIm(u) 2 -2v / 2MIm(ii)(Im(s)+Im(s / )) 



2lm(u)VMTT 



= 0, 



and, 



e -iMlm(uY e 2M(l-u z ) u AM 



_ e 2M(l-M 2 +21og|u|)_ 



Since \u\ < 1, we have 1 — \u\ 2 + 2 log \u\ < 0, and therefore 



lim e -2V2MIm(«)(Im(s)+Im( s ')) 
M— >oo 



e -4MIm(ti) 2 e 2M(l-« 2 ) u 4M 



We conclude that 



M 



lim DS 2M (uV2M + s, uV2M + s') = 0. 



And, since IS 2M (z,z') = -DS 2M (z, z'), 



lim IS 2 m(uV2M + s, uV2M + s') = 0. 

M^oo 



Turning to S 2M (uV2M + s, u\p2M + s'), we set 
Vm(s) = exp 

From Theorem [SI we have 



Vm{s)vm(s')S 2 m{uV2M + s, uV2M + s') 



2iV2Mlm(u) + {s- s') 



x exp (2iv / 2MIm(u)Re(s)) exp ( - 2iV2Mlm(u)Re(s')) 
x exp ^4Af(Im(w)) 2 - 2iV2MIm(u)(s - f ) - i(s - s' f 

^erfc |2\/MIm(w) + V2Im(s)| erfc \2\/MIm(u) + V2Im(s')| 
x exp (-2M |u| 2 - (su + s'u)\/2M - ss'^j e M (2M\u\ 2 + [su + s'u)^/2M + ss') 

Using (|9.12[) . we see 



1 



1 



Vm(s)vm(s')S 2 m(uV2M + s, uV2M + s') ~ - exp ( -~(s - s') 2 - M s ) 2 - Im(s') 2 

7T \ 2 

x exp (-2Af |u| 2 - (su + s'u)V2M - ss'^j e M {%M \u\ 2 + {su + s'u)V2M + ss') 
And thus, by Lemma [ 



1 



7T 
1 



1 



r/ M (s)VM(-s')S 2M (uV2M + s, uV2M + s') ~ - exp ( -^(s - s') 2 - Im(s) 2 - Im(s') 2 



exp 



M 2 W\ 2 
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Next, we set Dm to be the 2m x 2m diagonal matrix given by 

T> M = diag(n M (si),riM(-si), . . . ,r)M(s m ),r){-s m )) , 
noting that detDj^f = 1. It follows that 
Jim Ro tm (-,z) = lim Pf (D M [# 2 M(sk,s fc /)]!™ D M X 

M— too M— >oo \ ' — - 1 - 



= Pf 



= det 



_Iexp (-^ 



J -"V L + s kSk 



iexp (-^ - 



k,k' = l 



1 

— exp 



k/d 2 M 2 



+ s k s k 



k,k' = l 



□ 



where the last equation follows from Section 4.6 of [ST]. 

In order to prove Theorem [10] we also need to analyze the large M asymptotics of 
tm{u^/2M + s, u\[2M + r) where r and u are real numbers with < \u\ < 1 and s is a 
complex number. 

Lemma 9.3. Let r and u be real numbers with < |u| < 1 and let s be a complex number 
in the closed upper half plane. Then, 

lim r M (uV2M + s, uV2M + r) = 0. 

M->oo 



Proof. 

r M (uV2M + s,uV2M + r) = sgn(uV2M 



D -s 2 /2 2 M - 3 / 2 e~ M '" 2 



-usV2M 



fa r(2M - 1) 



x (uV2M + s) 2M_1 • 7 (m - ^,Mu 2 + uV2Mr + T — . 



We simplify this using Legendre's duplication formula for T(2M — 1) and by setting P(a, x) 
l{a,x)/T(a). 



-s 2 /2 

r M (uV2M + s, uV2M + r) = sgn(uV2M + ry 



V2 



-uV2Ms 



-7-)' 

uV2M J 



2M-1 



M M+l/2 e -Mu 2 U 2M-1 

Mi 



P( M - -^,Mu 2 + uV2Mr + — 



Next we use Stirling's approximation for M\ in the denominator in to show that 



r M (uV2M + s, uV2M + r) ~ sgn(u) 



-s 2 /2 



-uV2Ms 



2uy / TT 



exp (M(l - u 2 +21og|u|)) ■ P[ M - -,Mu 2 +uV2Mr 



uV2M J 
1 



2M-1 



Using the fact that 



2M-1 



sy/2M/u -s 2 /2u 2 
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wc find 



r M (uV2M + s,uV2M + r) 



sgn(^) 



ex P rr 



exp sV2M 



1-u' 



(9.13) 



/ 1 , r 2 

x exp (M (1 - u 2 + 2 log |u|)) • P ( M - - , Mu 2 + uV2Mr + — 

Finally, we notice that if < \u\ < 1 then 1 — u 2 + 2 log |u| < 0. It follows that 

lim exp u — u )) exp(M(l-u 2 + 21og|w|)) =0, 

and the lemma follows since 

< p(m - i,Mu 2 + uV2Mr + y^j < 1. 

Proof of Theorem \1(A First we consider the case where s and s' are both in the open upper 
half plane. From Theorem [U 



□ 



S 2 m(uV2M + s, uV2M 



+ s' 



2tt 



(s' - s)y erfc (V2Im(s)) erfc (\/2Im(s')) 



x exp ( - 2Mu 2 - uV2M (s + s') - ss')e M (2Mu 2 + uV2M (s + s') + ss 1 ) . 
DS 2 m{uV2M + s,uV2M + s') = ^= ) (s^ -s)^/erfc(v / 2Im(s))erfc(72Im(5 / )) 



/27T 



x exp ( - 2Mu 2 - uvWfs + s') - ss')e M (2Mu 2 + H\/2M(s + s') + ss'), 



and 



75 2 m(uV2M+s,uV2M + s') = 



_ e -L(s-s'Y 



-{s' ~ a) J erfc (v^Im(s)) erfc (-\/2Im(s')) 
x exp ( - 2Mu 2 - uV2M (s + s') - ss') e M (2M u 2 + W2M (s + s') + ss') . 

By Lemma 19.21 these converge to the appropriate entries of the complex/complex kernel 
K(s, s') as M -> oo. 

Next we turn to the case where u and r are real, and s is in the open upper half plane. 
In this case, Theorem [5] yields, 

~ ie~5( r - ¥ ) 2 / — 

S 2 m(uV2M + r, uV2M + s) = j= — (s - r)v/ erfc (V2Im(s)) 



/27T 

x exp ( - 2Mu 2 - uV2M(r + s) - rs)e M (2Mu 2 + uV2M(r + s) + rs), 

e -i(r-sf 



S 2 m(uV2M + s,uV2M - 



-r = 



12m 



erfc (V2Im(s)) 



x exp ( - 2Mu 2 - uV2M(r + s) - rs)e M {2Mu 2 + uV2M(r + s) + rs), 
+ r M {uV2M + r, uV2M + s). 
and thus, by Lemmas 19.21 and 10 . 31 



M 



lim S 2M (uV2M + r, uV2M + s) 



271 



(s -r)J erfc (V2lm(s)) , 
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and 



lim S 2 m(uV2M + s, uV2M + r) 

Af-»oo 



2tt 



erfc (V2lm(s)) . 



The limiting values for DS 2M (uV2M + r,uy/2M + s) and IS 2M (uV2M + r,uV2M + s) 
follow from this as well, since 

DS 2M (z,z') = -iS 2 M{z,^) and IS 2M (z,z') = iS 2M {z',z). 

Finally, we turn to the case where u, r and r' are all real. Here, Theorem |8] implies that 



S 2M (uV2M + r, uV2M + r') 



e 2 



i(r-r') 2 



x exp ( - 2Mu 2 - uV2M(r + r') - rr')e M (2Mu 2 + uV2M(r + r') + rr') 
+ r M (uV2M + r,uV2M + r'). 

and 

DS 2M (uV2M + r,uV2M + r') = ^ ~T^ e~g^~^' 

x exp ( - 2Mu 2 - uV2M(r + r') - rr')e M (2Mu 2 + uV2M(r + r') +rr'), 
From Lemmas 19.21 and 10.31 we see that 



1 



lim S 2M (uV2M + r, uV2M + r') = 

M->oo J2-K 



and 



lim DS 2M (uV2M + r, W2Af + r' 



= (^) e -| (r -o ! 
v/2^ 



All that remains to show is 



~ ; , l /Ir-r'l 

lim IS 2M (uV2M + r, uV2M + r') = - sgn(r' - r) erfc - — 7 =- A 

M— «x> 2 V v2 



First we write 



IS 2M (x,x') 



72 



2V? 

e -' 2 /2 



sgn(x') |/(a;,a;') - 
sgn(:r) \ I(x',x) 



*"/2 e -t 



-j=-C M {xV2t)dt 



2^ 



72 e -t 



C M (x'V2t)dt 



where I(x, x') is given as in (|9 . 1 0(1 and Cm = cosh— cm- That is, 

~ , e~ x2/2 e~ x ' 2/2 

IS 2M {x, x) = t= sgn(x')/(x, x 1 ) ^- sgn(x)/(x', x) 

2v 7r 2V7T 



• sgn(x) 



72 e -t 



-x 2 /2 

,_ C M (a/V2t) dt - sgn(>') 



*' 2 /2 e -t 



C M {xV2t)dt 



(9.14) 
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Making the substitutions x = u\j2M + r and x' = u\/2M + r', and assuming that M is 
sufficiently large, (|9.11|) yields 



sgn(^) 
2^ 



■ erfc 



I(x, x) = < 



4 

e x 2 /2 



V2 



■ erfc 



V2 



if u > 0; 



if u < 0. 



Consequently, 



lim 



-s 2 /2 e -^' 2 /2 

sgn(x')/(a;, x') =- sgn(x)I(x' , x) 



2^ 



- erfc 
4 



1 r f r — r' 
= — erf 
2 



erfc 

4 



V2 



V2 

Thus, if we can show that the second line of (|9. 14|) goes to as M — ► oo, we will have 

1 „ I ' t — r n 



M 



as desired. 

We thus consider 



Clearly, for any v > 0, 



lim j - sgn(x - x') + IS 2 m(x, x') \ = ^ sgn(r - /) - ^ erf 

f^oo Z J Z Z 



= — sgn(r — r') erfc 



V2 
\r-r'\ 

~7T~ 



2^ 



sgn(x) j -^C M {x'V2t)dt 



C m (v) < — = — - - ^ 1 • 



x m! T(2M) 



Thus, 



It follows that 



* 2 /2 P -t 



^C M (x'V2~t) dt< —dt = V^erf [ J=L < V^. 



v/2 



e -x' 2 /2 rx 2 /2 - t e~ X ' 2 / 2 

and making the substitution x' = u\/2M + r', this goes to as M — ► oo. 



□ 



9.5 The Proofs of Theorem [12] and Theorem 



In order to prove Theorem [12] we need analogs of Lemmas 19.21 and 19.31 for the case where 
u 2 = 1. 
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Ginibre ensemble of random matrices and its scaling limits 



Lemma 9.4. Suppose a is a complex number and 

, M = 1 + 7 J= + (M-), 

Then, 

lim e- 2Mv -e M (2Mv M ) = ~ erf c ( . 
Proof. As in the proof of Lemma 19.21 we set v = Vm and write 

w = WM = v + W—v 

so that 

e- 2Mv e M (2Mv) = e~ {2M ~ 2) ^ 'e M ((2M - 2) W ). 

We are now in position to use a result of Bleher and Mallison [3l Theorem B.l], which 
shows that 

e -(2M-2)«, eM ( (2M _ i)w) „ 1 erfc ^( w)v / (2 M-2)) , (9.15) 

where 

{(w)= fu^_(^ (u^ _ 
\/2 6\/2 36V2 

In our case, 

eH = -^= + 0(M- 1 ), (9.16) 



and the Lemma now follows from (I9.15P and (|9.16j) . □ 

Lemma 9.5. Let u = ±1, r€R and let s be in the closed upper half plane. Then, 

1 _ 2 

lim Tm (u\ 2M + s, uy 2M + r) = — i=e s erfc(— ur). 

M->oO 4^/7T 

Proof. From (|9.13[) we have that 

, , 1 2 / 1 , r 2\ 

r M (uV2M + s,uV2M + r) = ^-^e' 8 P [M — -,M + uV2Mr + — (9.17) 



2y^ 

In [25], Temme gives the uniform asymptotic expansion for P(a,x) when a > and x £ 

P{a, x) ~ - erfc (sgn(l - X)y/ a(X - 1 - log A)) ; A = — . 

2 V /a 

In our situation, 

A= l+^ + ^ =1 + ^ + 0(M _ 1}- 



1 - VM 



It follows that, as M — » oo, sgn(l — A) — » — usgn(r), and 



Thus, 



A-l-logA= — + 0(M"^). 



I , r 'A i 

P( M , M + uV2Mr + — J ~ - erfc (-ur) . 



and the lemma follows from (|9.17[) . □ 
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Proof of Theorem WA The proof of Theorem [12] is the same, mutatis mutandis, as that of 
Theorem [TU] replacing the asymptotics in Lemmas 19.21 and 19.31 with those in Lemmas 19.41 
and EH □ 

Proof of Theorem \1S\ The proof of Theorem [13] is the same, mutatis mutandis, as that of 
Theorem QT] replacing the asymptotics in Lemmas 19.21 with those in Lemmas 19.41 □ 



Appendices 



A Correlation Functions for (3=1 and (3 = 4 Hermit ian 
Ensembles 

In this appendix we will use the Pfaffian Cauchy-Binet Formula (see [B]) in order to derive 
the correlation functions of the [3=1 and (3 = 4 Hermitian ensembles. We will keep 
the exposition brief, but will introduce all notation necessary for this appendix to be read 
independently from the main body of the paper. We reuse much of the notation from main 
body of the paper so that we may also reuse the same proofs. For convenience, N will be a 
fixed even integer; similar results are true for odd integers. 

Given a Borel measure v on R we define the associated partition function to be 



^ ■ JR« 



where A('y) is the Vandermonde determinant in the variables 71,72, ■ ■ ■ ,7iv and vm is the 
product measure of v on M. N . When (3 = 1 we define the function £ : R 2 — > {—^,0, i} and 
the operator e v on L 2 (u) by 

£(7,7') := -sgn(7-7') and e v g{i) := / g(y)£(y,l) du{y). 

When (3 = 4 we define £(7,7') := and e v g(y) ■= g'(y)- We use e„ to define the skew- 
symmetric bilinear form on L 2 (u) given by 

(9\h) v := f (<?(7)^(7) - «*ff(7W7))<M7). 

JR 

Theorem A.l. Let b := \ffi, and let q be a family of Nb monic polynomials such that 
deg q„ = n . Then, 

where Uq = [{q n \q n >)v\\ n, n' = 0, 1, . . . , Nb - 1. 

This theorem follows from de Bruijn's identities [7]. 

We set A to be Lebesgue measure on R. If there is some Borel measurable function 
w : R — * [0, 00) so that v = wX (that is, dv/d\ = w) then we define Z w := Z v . Clearly, 

ZW = Jn.J RN nNil)dXNh) Wh6re ^(7):={lI^(7n)}|A(7)l /3 - 
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We may specify an ensemble of Hermitian matrices by demanding that its joint probability 
density function is given by jv - The nth correlation function of this ensemble is then 
defined to be R n :!"-> [0, oo) where 



R n {y) 



i 



i 



(N - n)\ 



^jv(y V 7) d\ N - n {~f), 



(A.l) 



where y V 7 € is the vector formed by concatenating the vectors y € W 1 and 7 € R" - ". 
By definition, Rq = 1. Here we take (|A.1[) as the definition of the nth correlation function; 
one can use the point process formalism to show that this definition is consistent with the 
definition derived in that manner. See [4] for details. 

We set /Lt n)n ' to be the n,n' entry of (Uq A ) _T , and we define q n := wq n . Using this 
notation we define the functions Sn,ISn and DSn 

Nb-l 

b 



by 



Nb-l 



2 ^ 2 . 

SW(7,7') : = T ^2 Vn,n> qn(l)e\q n >(l'), IS N (l , l') ■= T ^ Vn,n> t\q n {l) t\q n > (l' ) 



Nb-l 



and 

2 , 

DS N {~f,~f') \= - ^ Mn,n'5n(7)?n'(7') 
n,n'— 

The matrix kernel of our ensemble is then defined to be 



-S N W,<y) 1^(7,7') + £(7, 7')- 



Theorem A. 2. 



Rn (y) = Pf [K N ( Yj , yj,)] ; j, j' = 1, 2, . . . , 71. 

Our proof of this theorem begins by setting r\ to be the measure on M given by 

T 



where 2/1 , 2/2 , ■ • • , 2/t are real numbers and Ci, C2, • • • , or are indeterminants and 8 is the prob- 
ability measure on K with point mass at 0. We will assume that T > N. As with Theorem [4] 
in the main body of this paper, we will prove Theorem I A. 2 1 by expanding Z w ^ x+ri " > / Z w in 
two different ways and then equating the coefficients of certain products of Ci, c%, . . . , cy. 



Proposition A. 3. 



= Pf (J + [V^K N (y t , y t ,)]); t, t' = 1, 2, . . . , T, 



where J is defined to be the 2T x 2T matrix consisting of 2 x 2 blocks given by 



J : = 





' 1" 




S t ,t> 


-1 





M = l,2,. 



Proposition lA.3l is proved in Section [9Tl 

For each N > we define 3^ to be the set of increasing functions from {1,2, ... ,n} into 
{1,2, ... ,T}. Given a vector y g R T and an element t £ 3^, we define the vector yt £ K ra 

by y t = {jA(i),2/t(2), • ■ -,yt( n )}- 
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Proposition A. 4. 



N , T . 

= 1 +EE n^K^)' 

n=l ie3 T ^ .7 = 1 J 



(A.2) 



tea? 

The proof of Proposition IA.4I is given in Section 19.21 
Finally, we set K to be the 2T x 2T block matrix given by 

K := y^K N {y u y v )]] t, t' = 1, 2, . . . , T. 

From the formula for the Pfaffian of the sum of two antisymmetric matrices (see Proposi- 
tion \7§ and Proposition IA. 31 we have that 

T 

Pf[J + K] = l + £ E PfKt - 

n=i teaj- 

where for each t G 3%, K t is the 2n x 2n antisymmetric matrix given by 



Finally, 



K t = [y/cuj) Ci( f )K N (y i{j) , y l{r) )] ; j,f = 1,2,...,S. 
Pf K t = { f[ c t(i) | Pi[K N (y m ,y l{jl) )]; j, f = 1,2,... 

L 3=1 J 



(A.3) 



and Theorem IA.2I follows from Proposition IA.41 by comparing coefficients of c\c 2 ■ ■ ■ c n in 
(|AT2|l and CO)) . 



B The Pfaffian Cauchy-Binet Formula 

Theorem B.l (Rains). Suppose B and C are respectively 2J x 2J arid 2i4T x 2K antisym- 
metric matrices with non-zero Pfaffians. Then, given any 23 x 2K matrix A, 

Pf (C- T - A T BA) _ Pf (B- T - ACA T ) 
Pf(C- T ) ~~ Pf(B" T ) ' 

Proof. Let I 2 k and I 2 j be respectively the 2K x 2K and 2 J x 2 J identity matrices, and let 
O be the 2 J x 2K matrix whose entries are all 0. Then, an easy calculation shows that 



" hj 


O " 




B T 


A" 




'I 2 j 


BA" 




A T B 


^-2K_ 




A T 


c- T _ 




o T 







B~ T 

o T 



o 

C T A T BA 



and similarly 



~h.I 


-AC" 




B T A" 




hj o " 




B" T -ACA T O 


o T 


l2iY 




A T C- T _ 




CA T I 2K _ 




o T C T 



Now, if D and E are 2N x 2N matrices and D is antisymmetric, then it is well known that 
Pf (EDE T ) = Pf D • dctE. From which we conclude that 



Pf 



B T 

o T 



O 

C T A T BA 



= Pf 



B T ACA T 

o T 



O 



The theorem follows since the Pfaffian of the direct sum of two even rank antisymmetric 
matrices is the product of the Pfaffians of the two matrices. That is 

Pf(B" T ) Pf(CT T - A T BA) = Pf(B" T - ACA T ) Pf(CT T ). □ 
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C Limiting Correlation Functions for the Complex Gini- 
bre Ensemble 

The complex Ginibre ensemble consists of N x N complex matrices with i.i.d. normal entries. 
In this section we derive the scaling limits of the correlation functions of the Ginibre complex 
ensemble in the bulk and at the edge. As the quantities of interest are similar to those in 
the main body of the paper we will reuse much of our previous notation for the analogous 
quantities. 

In his original paper on the subject, [13], Ginibre showed that the joint density of eigen- 
values is given by 

^ ( N 1 
W (7) = ^ Rwhn) |A(7)| 2 , 

where w(j) — e™' 7 ' , A (7) is the Vandermonde determinant whose columns are given in 
terms of 71, 72, ... , 7jv, and 

We may take the nth correlation function of this ensemble to be the function R n : C" — > 
[0, 00) given by 

1 1 f 

Z (N — n)\ JcN-n 

The correlation functions can also be defined as densities with respect to Lebesgue measure 
which satisfy an identity analogous to ()5. 1|) . 

Ginibre gave a closed form for R n in terms of a scalar kernel. Specifically, he showed 
that 

R„{z) = det [ifiv(zfc,2fe')]fc,fc'=i ' 

where 

K N {z,z>) = i-cxp (-^ - ¥p) e M {zn 

Clearly then, 

lim R n (z) = det 

This is the limiting correlation function of the complex Ginibre ensemble at the origin. 
Notice that an almost identical expression appears in Theorem I 111 

Like the real Ginibre ensemble, the complex Ginibre ensemble satisfies the circular law. 
We therefore expect that limiting correlation functions will emerge after scaling eigenvalues 
by a factor of ^/N. 

Theorem C.l. Let u be in the closed unit disk, and suppose s\, S2, • ■ ■ , s n are complex 
numbers. Set 

z k = uvN + Sk k = 1, 2, . . . , n. 



1 / \z k \ 2 W>\ 1 M _ \ ' 



Then: 
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1. Limiting correlation functions in the bulk. If \u\ < 1, 

1 



lim R n (z) = det 



2tt 



exp 



\s k \ 2 M 2 



SkSk> 



fc,fc' = l 



2. Limiting correlation functions at the edge. If \u\ = 1 

1 



lim R n (z) = det 



2tt 



exp 



s fc | 2 |s fc '| 2 _ \ (s k u + s k ,u 
— h SfcSfc' erfc 



k,k' = l 



Proof. Let AT = 2M , 



, (su — su) /— 
i) M (s) = exp ( v ^_ > VM 



and define D to be the n x n matrix given by D = diag (i/>m(si), (§2), ■ ■ • , ^(s ti )). It is 
easily seen that \4>m(s)\ = 1 and detD = 1. We also define K to be the n x n matrix given 

by ' 



K = 



1 

2^ 



exp 



Sfe' 



s k s k > erfc 



SkU + Sfc'U 



k,k' = l 



Then, 



lim R n U) = lim det (DKXT 1 ) 



(C.l) 



The k, k! entry of K can be computed to be 



— exp -2M\u 
2tt 



2 (s k u + s k 'u) {s k u + Sk'U) |s fc | 2 \sk'\ 2 



V2 



V2 



e M (2M\u\ 2 + {s k u + s k >u)V2M + s k s k ^ 



1 ( o«fi 12 1 -.- \ /^77 - ^ f l^l 2 l^'l 2 , _ 
= — exp \ -2M\u\ - (s k u + s k 'u)V2M - s k s k ' j exp I h SfcSfc< 

f{s k U-S k u) j—\ ( -(Sk'U-Sk'U) r— 

x exp ( ^ VMj exp Vm 

x e M (2M|m| 2 + (s k u + s k >u)V2M + s k s k >^J . 

It follows that the k, k' entry of DKD 1 is given by 
1 



2tt 



exp 



\s k \ 2 \s k >\ 2 



SfeSfe' exp -2M\u\ - (s k u + s k >u)V2M - s k s k 



e M (2M\u\ 2 + (s k u + s k 'u)V2M + s k s k ^ 



Statement Q] of the theorem now follows from (|C.1[) and Lemma 19.21 Statement [2] follows 
from (jC.lj) and Lemma l9~4l □ 
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D Plots of Correlation Functions for the Real Ginibre 
Ensemble 

This appendix contains various visualizations of the limiting correlation functions of the real 
Ginibre ensemble. 

As usual, H (respectively H) is the open (closed) upper half plane. Given a point u € H 
with \u\ < 1, n, 7*2, . . . , rt € R and s\, S2, . . . ,s m € H m , we set 

= u V 2M + rj j = 1,2, and Zk = uV2M + Sk k = 1,2, . . . ,m. 

We will use the notation 

H% m (ri,...,rt,si,...,s m ) = lira i? f , m (x,z), 

M— >oo 

where i?^ m (x, z) is the £, m correlation function of the real Ginibre ensemble of 2M x 2M 
matrices. 

D.l The real bulk 

Let u be a point in the real bulk. The local density of real eigenvalues is constant (equal 
to 1/^/2%). The limiting correlation function, R% (ri, r-i) is invariant under real shifts, and 
hence i?2o( r ii r 2) = -^2o( r i — r 2>0). We may therefore plot this correlation function as a 
function of a r\ — r%. As |ri — ra| — > oo, this quantity approaches (2tt)^ 1 — the square of the 
density of real eigenvalues. See Figure [T] 

The local density of complex eigenvalues in the real bulk is given by Rq 1 (s). Due to the 
invariance of the correlations functions with respect to real shifts, this is a function of Im(s) 
only. As Im(s) — * oo this density tends toward the density of eigenvalues in the complex 
bulk. Specifically, 

km R% A (s) = -. 

Im(s)— >oo 7T 

See Figured 

The correlation function i?" x(r, s) is invariant under real shifts, and thus can be plotted 
as a function of r — Re (s) andlm(s). As Im(s) — > oo, R^^.s) approaches 2 -1 / 2 7r~ 3 / 2 — the 
product of the density of real eigenvalues in the real bulk and the density of eigenvalues in 
the complex bulk. See Figure [31 

D.2 The complex bulk 

When u is in the complex bulk; that is u E H and \u\ < 1, the density of eigenvalues 
is constant (equal to 1/tt). The only non-trivial correlation function we can visualize is 
Rq 2 (s, s'). This function is invariant under both real and imaginary shifts. That is, we may 
plot Rq 2 (s, s') as a function of s — s' . As \s — s'\ — > oo, -Ro.2( s ; s ') approaches 1/ir 2 — the 
square of the density of eigenvalues in the complex bulk. See Figure [4] 

D.3 The real edge 

For concreteness we will concentrate on the real edge corresponding to u — 1. At the real 
edge the local density of real eigenvalues is no longer constant. Here the density is given 
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by R\ o(r). As r — > — oo we expect the local density of eigenvalues to approach the density 
of real eigenvalues in the real bulk, l/y/2n. Indeed, this is the case. As r — > oo the local 
density of eigenvalues decreases to 0. See Figure [5] 

The local density of complex eigenvalues at the real edge is given by Rq 1 (s). This can be 
plotted as a function of Re(s) and Im(s). If s 'moves' in the direction of the complex bulk 
(loosely speaking, Re(s) — > — oo while simultaneously, Im(s) — > oo) then Rq 1 (s) approaches 
1/tt — the density of eigenvalues in the complex bulk. Since the real axis repels complex 
roots, Rq 1 (s) approaches with Im(s). See Figure 

We may also plot R\{r, r'). Here we expect that if r — * — oo and r' — ► oo and |r— r'| — ► oo, 
then R\(r, r') should approach the square of the density of real eigenvalues in the real bulk, 
(27r) _1 . If r — r' approaches 0, then the repulsion of eigenvalues implies that R\{r, r') — > 0. 
Similarly, if either r or r' is large and positive, then we are looking at the local correlation of 
involving a real eigenvalue away from the bulk and therefore R\{ r i r ') is small. See Figure[7] 

D.4 The complex edge 

The limiting kernel at the complex edge, when u is on the unit circle, is invariant under 
shifts in the direction of the tangent line of the unit circle at u. For simplicity we take u = i, 
so that, for instance, the local density of complex eigenvalues is invariant under real shifts. 
This density if given by Rq 1 (s) which is a function of Im(s) only. As Im(s) — > oo we are 
moving away from the bulk and thus R l ^(s) — * 0. If Im(s) — ► — oo then R 1 (s) approaches 
2/tt, the density of complex eigenvalues in the bulk. See Figure [8] 
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-x.32 




Im(s) 



Figure 2: The density of complex eigenvalues in the real bulk as a function of Im(s). 



Figure 3: R\ t i(r, s) in the real bulk as a function of r — Re(s) and Im(s) 
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Figure 6: The density of complex eigenvalues at the real edge as a function of Re(s) and 
Im(s). 



4G 




Figure 8: The radial density of eigenvalues (represented here as R l Q1 (s)) as a function of 
the radius (given here by Im(s)) at the complex edge. 



Alexei Borodin 

California Institute of Technology 

Pasadena, California 

email: borodin@caltech.edu 

Christopher D. Sinclair 
University of Colorado 
Boulder, Colorado 

email: chris . sinclairOcolorado . edu 



